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The OZI rule is prominent in hadronic phenomena only because OZI vio- 
lation is typically an order of magnitude smaller than expected from large Nc 
arguments. With its standard ^Pq pair creation operator for hadronic decays 
by flux tube breaking, the quark model respects the OZI rule at tree level 
and exhibits the cancellations between OZI-violating meson loop diagrams 
required for this dramatic suppression. However, if the quark model explana- 
tion for these cancellations is correct, then OZI violation would be expected 
to be large in the nonet with the same quantum numbers as the pair cre- 
ation operator: the 0'*''^ mesons. Experiment is currently unable to identify 
these mesons, but we report here on a lattice QCD calculation which confirms 
that the OZI rule arises from QCD in the vector and axial vector mesons as 
observed, and finds a large violation of the rule in the scalar mesons as antic- 
ipated by the quark model. In view of this result, we make some remarks on 
possible connections between the ^Pq pair creation model, scalar mesons, and 
the anomaly responsible for the large OZI violation which drives the r/' 

mass. In particular, we note that our result favors the large Nc and not the 
instanton interpretation of the solution to the r]' mass problem. 
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I. BACKGROUND 



The phenomena which led to the formulation of the OZI rule have had a definitive 
impact on our understanding of strong interactions. The fact that "aces" (i.e., quarks) led 
to a simple interpretation of the properties of the meson was clearly a very important clue 
for Zweig since it was natural for the (f) to be pure ss and for certain production cross 
sections to be small so long as "hairpin graphs" were dynamically suppressed (see Fig. 0). 

The dynamics behind the suppression of hairpin graphs in QCD has remained unex- 
plained. The phenomenology of meson mixing angles in QCD-based quark models was 
described in the mid-1970's in a number of papers In such models, processes with 

the quark line topology of the double-hairpin graphs of Fig. ^b) (but with arbitrary time 
orderings) modify the quark- ant iquark transition amplitudes from the totally flavor diagonal 
form associated with the "scattering" quark line topology of Fig. ||(a), namely 
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(for illustrative purposes we have suppressed all space-time labels and specialized to the 
case of SU (2) flavor where the matrix spans the basis ud, du, uu, dd) by the addition of the 
annihilation amplitudes A 
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Using this framework 0, it was noted that the OZI mixing amplitude A characterizing Fig. 
H(b) was of order 10 MeV in the established meson nonets, with the sole exception of the 
ground state pseudoscalar meson nonet, where A is an order of magnitude larger. These 
observations were consistent with the pattern one would expect for heavy quarkonia where 
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the ground state pseudoscalar double-hairpin is larger than the vector double-hairpin by one 
factor of (tts/Tr)"^, and excited state double-hairpins are suppressed by having vanishing 
wave functions at r = 0. However, an explanation for this pattern in light quark systems 
was lacking. 



FIG. 1. A typical hairpin reaction, where the two lines-ellipsis-two lines on the left denote an 
arbitrary OZI allowed process and an ss hairpin is shown for concreteness. Note that gluonic fields 
and closed qq loops are not represented since the external quark line topology is all that is relevant 
to the rule. 

The large size of the ground state pseudoscalar double-hairpin is a manifestation of 
the "f/A(l) problem" 0: the equations of motion of QCD, taken naively, would imply that 
spontaneous chiral symmetry breaking leads to nine and not just eight Goldstone bosons P], 
but the large mass of the rj' seems to disqualify it from the role of the flavor singlet Goldstone 
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boson. However, the Ua{1) current is anomalous, and by the late 1970's it was understood 
through the study of instantons P,p!0[] that the anomaly leads to a nonconservation of the 
Ua{1) charge and thereby to the evasion of Goldstone's theorem in the flavor singlet channel 
when chiral symmetry is spontaneously broken. The connection between the quark model 
picture of double-hairpins and instantons was discussed by Witten lO, Veneziano [IT2 



and others, who explored more generally the conflict between instantons and the large Nc 
expansion [ |T^ . 
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(a) (b) 

FIG. 2. Quark line diagrams associated with the OZI rule in meson nonets: (a) a normal 
OZI-conserving quark-antiquark "scattering" process, (b) a typical double-hairpin "annihilation" 
process leading to flavor-mixing in meson wave functions. Since the OZI-violating reactions of Fig. 
1^ could occur through OZI-allowed meson emission followed by flavor mixing, these processes are 
the simplest manifestations of OZI violation. Note that gluonic fields and closed qq loops are not 
represented since the external quark line topology is all that is relevant to the rule. 

4 



(The reader familiar with instanton lore may be puzzled by the connection between 
the annihilation amplitudes Aq^j of Eq. (2) in the pseudoscalar mesons and instanton- 
induced effects in the pseudoscalar mesons. The latter effects are associated with the 't 
Hooft interaction [|T^ which (in our illustrative SU{2) flavor case) leads to uu —>■ dd and 
dd uu but not the diagonal entries in Eq. (2) for AT corresponding to uu uu or 
dd dd transitions. Recall, however, that the 't Hooft interaction also has ud ud and 
du — > du interactions, i.e., the S'-like amplitudes of Eq. (1). Thus the instanton-induced 
interactions also admit the decomposition of Eqs. (1) and (2) with S = —A. We will 
elaborate upon this point below.) 



TABLE I. OZI-violating amplitudes in meson nonets. These amplitudes are defined to be the 
contribution of the uu dd double-hairpin to the nonet mass matrix. 



nonet 


empirical 
(MeV) 


quark model loop 
contribution to Aq'^j (MeV)^ 


0-+ 


+400 ± 200 * 




1" 


+7± 1 


-2 + 4 


2++ 


-22 ±3 


+6 + 14 


1++ 


+11 ± 15 


+12 + 12 


0++ 


see text 


-450 + 200 * 


1+- 


-32 ± 12 


-15 + 7 


3~ 


-12 + 4 


+4 + 7 


4++ 


+6 + 18 


+16 + 7 



See Ref. [HI. 



^The quoted "theoretical error" assigned here is the range quoted in Ref. [|T^] for meson loop 
processes from reasonable parameter variations. 
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The large Nc expansion is the only known field-theoretic basis for the general success 
of the valence quark model, Regge phenomenology, the observed narrowness of resonances, 
and the OZI rule. In particular, the OZI-violating meson mixing amplitudes of Fig. ^(b) 
are all of order l/Nc- Ironically, such a suppression of these amplitudes seems perfectly 
consistent with the effects in the pseudoscalar mesons, but not strong enough to account for 
the extremely small amplitudes seen in other nonets. See the second column of Table |. 

The unexpected suppression of most OZI-violating amplitudes beyond a simple factor 
of 1/Nc is elevated from a dynamical puzzle to a paradox when the various time-orderings 
of Fig. H(b) are projected into a hadronic basis. In such a basis, fiavor mixing could arise 
through an intermediate glueball, through an instantaneous interaction, or via a hadronic 
loop process in which Fig. ^(b) has the time-ordering shown in Fig. |^. The paradox arises 



from the observation that these OZI-violating hadronic loop processes can proceed by 
sequential OZI-allowed vertices with known and unsuppressed strengths. These hadronic 
loop diagrams may be associated with contributions to meson propagators arising from 
second order (real and virtual) decay processes, and as such are of order (l/y^iVc)^, as 
expected. This factor of 1/Nc is also perfectly consistent with the observation that the 
imaginary parts of these propagators give the l/Ai'c-suppressed meson widths which are 
generally of order of hundreds of MeV. Nevertheless, OZI phenomenology requires that the 
1 /A^c-suppressed real parts (from the full meson spectrum and not just the kinematically 
allowed part) be an order of magnitude smaller. Explicit model calculations substantiate the 
generic result that individual hadronic channels of the type depicted in Fig. ^ would indeed 
contribute hundreds of MeV to OZI-violating meson mixing. Thus even if the other possible 
sources of OZI violation (from the other time-orderings) were dynamically suppressed, these 
hadronic loop diagrams would seem to spoil the OZI rule. This rule requires that Aqzi << 
— SO that the nonet and not the SU (3) limit is realized. Thus it is a necessary (but 
not sufficient) condition for the OZI rule that there be some conspiracy between hadronic 
loop processes which suppresses them below their expected 1/Nc strength |TH . 
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FIG. 3. Two sequential OZI-allowed processes can lead to the topology of Fig. |2|(b). 



II. A PROPOSED RESOLUTION 



The authors of Ref. [|T7| proposed a resolution of this paradox. They examined the OZI- 



violating amphtudes Aqzi in non-pseudoscalar channels from the complete tower of hadronic 
loop processes to determine if "miraculous" cancellations between the hundred-MeV-scale 
real parts of individual channels could be responsible for the suppression of the sum over 
channels beyond a simple power of l/N^. To make such a calculation one must have a 
complete model for meson trilinear vertices since, if such a conspiracy is to occur, it will have 
to be based on an underlying pattern of coupling strengths and signs. The nonrelativistic 
quark model is complete in this sense: using the standard ^Pq pair creation operator for 



hadronic decays by flux tube breaking |18| and valence quark model wave functions, all 



trilinear vertices and their associated form factors are prescribed. Since the model is a 
nonrelativistic one, only the time-ordering of Fig. R| with a two meson intermediate state 



can be calculated in this way, but within this framework Ref. [|T7| shows that in general 
a "miraculous" cancellation between channels does indeed occur. This cancellation occurs 
between groups of intermediate meson states that might have been difficult to anticipate 
a priori. Consider the prototypical case of u — (p mixing where u — {AB)l with 
L the AB relative angular momentum. The intermediate states contributing to Fig. § are 
{KK)p, {KK*)p, {K*K*)p, {K*Kl)s, {KKa,)s, {KKa,)D, {K*K,,)s, {K*K,,)d, {KKt,)s, 
{KK,,)d, {K*K,Js, {K*KtjD, {KK*)d, {K*K*)s, {K*K*)s, ... where K and K* are the 
ground state {£ = 0) pseudoscalar and vector mesons, and Kq, Ka^, K^,^, and are the 
first excited state {^ = 1) strange mesons with = 0^, l"*", l"*", and 2'^ which would be 
associated with the oq, cti, &i, and a2 octets in the SU{3) limit. (Note that the ellipsis 
denotes more highly excited intermediate states, including ones in which each leg of the 
intermediate state is excited, and that charge conjugate intermediate states are implied.) 
As expected on the basis of the previously described arguments, a typical channel in this 
sum contributes of order 100 MeV to Aq^j. However, intermediate states with the same 
total orbital angular momentum but opposite values of (—1)^ tend to cancel. Thus, for 
example, the {£a = 0,£b = 0)p channels with Ltotai = (-a + (-b + L = 1 all have the same 
sign, but they strongly cancel against the = 0,£b = 1)s+{(-a = '^Ab = 0)5 channels! 

The calculation is formidable. With standard quark model parameters the form factors 
are quite hard and complete convergence is achieved only after summing of order 10 thousand 
channels, corresponding to Ltotai — 10. With reasonable variations of standard parameters 
the contribution of an individual channel waxes and wanes, as does the speed of convergence. 
However, the underlying mechanism of the cancellation is simple and very robust: Aq^j is 
much smaller than its component pieces because of an approximate "spectator plus closure 
limit". This limit is illustrated in Fig. ^ which shows the standard ^Pq operator with 
J^*^ = 0^+ trying to create and then annihilate quark-antiquark pairs with J^'-^ = 1 . 
If a single two meson intermediate state is inserted into this diagram, it will project out 
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pieces of this amplitude of order l/N^as expected, but if the original (final) qq pair does not 
distort the J^^ of the produced (annihilated) pair (the spectator approximation) a complete 
set of intermediate states with a common energy denominator (the closure approximation) 



will give zero amplitude. Ref. shows that deviations from this "spectator plus closure 
limit" are naturally small, leading to the observed order of magnitude suppression of the 
loop contribution to A]^^^ relative to expectations. See Table |. The interested reader 



is referred to Ref. [T^] for a detailed explanation of the resiliency of this limit. This quark 
model solution to the "second order paradox" associated with the OZI rule also appears to 
justify the conspiracies between Regge trajectories required to explain the suppression of 



cross sections requiring "exotic" exchanges (e.g., those with isospin 2) [jl9[. Since "exotic" 
exchanges can occur by double Regge exchanges (analogous to the second order loop pro- 
cesses), only a conspiracy between exchanges (analogous to the conspiracy between loops) 
can give the observed suppression of such cross sections. 

While the order of magnitude suppression of the loop contribution to A^^/ is robust, the 
contribution of individual channels and the residue after the cancellations have occurred is 
model sensitive, so a prediction for the actual value of this amplitude cannot be made. This 
is not a great loss, however, since the accuracy of the model is very suspect: its dynamics is 
nonrelativistic, and it has ignored the Z-graph time orderings of Fig. |^. More significantly, 
any such amplitude would need to be added to the unknown pure glue and instantaneous 
contributions to the qq — * q'q' transition before being compared to experiment. Thus the 



important conclusion of Ref. |jT^ is the qualitative one that the "second order paradox" can 
be evaded. 

Ref. [0 confirms that Aqzi from meson loop diagrams is small in not only the vector 
mesons but in all other well-established nonets: those with J^*-^ = 2++, 1^^,3 , and 
4"'"^. The key, of course, is that the nonet J^'" must differ from that of the ^Pq pair creation 
operator |2y]. From this simple requirement follows a rather spectacular prediction: OZI 
violation should be very strong in the scalar meson nonet. 

The scalar mesons, and especially the isoscalar scalar mesons which would display the 
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effects of OZI violation, have been notoriously difficult to understand experimentally. Over 
the last thirty years the mass of the hghtest isoscalar scalar meson quoted by the Particle 
Data Group has varied between 400 and 1400 MeV, while the quoted width has varied 
between and 100 and 1000 MeV. (We have removed the /o(980) from this compilation under 
the presumption that it is a KK molecule, or this spread of values would be even wider.) 
The experimental status of the scalar meson nonet becomes even more obscure when one 
recalls that the hghtest glueball is expected to have J^'-^ — 0"^"^ and a mass around 1.5 GeV. 
One can only say with confidence that the experimental situation does not exclude that 
A^ozi is large. 

Fortunately, there is an alternative to checking this prediction of the quark model mech- 
anism against experiment. We can check it against calculations from lattice QCD. 
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FIG. 4. A graphical representation of the "spectator plus closure limit": in this limit the 0"*"+ 
pair creation operator cannot destroy or create a 1 pair. 
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III. OZI ON THE LATTICE 



A. OZI Violation in the Quenched Approximation 

Matrix elements of the type {0\T[q' (y)T'^^^ q' (y) q{x)T'^^^ q(x)]\0) , where r"^^*^ carries 
space-time indices which determine the J^*^ of the propagator being studied and q' 7^ q, 
describe OZI violation in mesons. In leading order in 1/Nc, such processes can proceed 
through diagrams of the type depicted in Fig. ^(b) (of which Fig. | is the particular time- 
ordering relevant to the hadronic loop diagrams), i.e., they receive leading contributions 
in the quenched approximation in which internal quark-antiquark loops are ignored. (Of 
course the accuracy of the quenched approximation can be questioned, but this is irrelevant 
to the main points of this paper, including checking a prediction of the quark model in which 
internal quark loops are also neglected.) 

In the absence of OZI violation, the cj-like -^{uu+dd) and 0-like ss sectors are segregated 
and each develops its own tower of meson excited states of each allowed J^'-^ . If the OZI- 
violating amplitudes Aq'^j in that channel are small, then in leading order they simply 

tPC 

shift the masses of each state by ~ ^oz/ create u — 0-like mixing with a mixing angle 

7-PC 

~ Aqzj/ Am where Am is the unperturbed mass difference between the uj- and 0-like states 

7 PC 

being mixed. In such circumstances the empirical value of Aq^j may be extracted from either 
the u — p-like mass difference or the u — 0-like mixing angle and compared directly with 
the quenched lattice amplitudes since the latter may be construed as correctly representing 

7 PC 

OZI-violation in the quenched approximation in lowest order perturbation theory in Aq^j. 

7 PC 

If Aqzj is strong, as in the pseudoscalar channel, the situation is more complicated. In 
such circumstances two new effects come into play: the masses of u- and 0-like states can be 
shifted strongly, so that their mixing angle may not be determined by their unperturbed mass 
difference, and treating the mixed propagator from qq q'q' in lowest order in Aq'^j may not 
be valid. The former effect is straightforward, but the latter can be complex. For example, a 

7 PC 

higher order treatment of Aq^j appears to be inconsistent with the quenched approximation. 
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as shown in Fig. |5|. However, the process depicted in Fig. |^ is one of a series of processes 
with internal quark loops which arise from repeated iteration of the quenched amplitude. 
Their effect and that of the diagonal mass shifts is to create a propagator matrix with entries 
corresponding to the quenched approximation; when diagonalized perturbatively this matrix 
gives the masses and mixing angles for weak OZI violation, but for strong OZI violation it 

7 PC 

may be diagonalized exactly, thereby summing the series of sequential applications of ^oz/- 

7 PC 

Another closely related possible complication is that a large Aq^j can create strong mixing 
with the glueball sector, requiring that the propagator matrix be enlarged yet further. 




FIG. 5. A contribution to OZI-violating meson mixing which is superficially inconsistent with 
the quenched approximation. 
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For the pseudoscalar mesons, the preceeding discussion of the effects of a large Aqzi are 

particularly significant. In the chiral limit with Aq^j = 0, the Ua{^) meson the rj' 

is also massless. As a result, the quenched OZI- violating amplitude of Fig. | will give Aq^j- 
sandwiched between two massless propagators, i.e., it will give an rj' propagator that looks 
nothing like that of a massive, S'f/(3)-flavor-mixed t]'. In this case, to even qualitatively 
relate the quenched amplitudes to nature one must extract Aq^j from them and add these 
amplitudes to the propagator matrix (the broken SU{3) analog of Eq. (1)) which one 
diagonalizes exactly. The resulting full propagator will have a massive S'[/(3)-flavor-mixed 
7]' which sums the single particle effects of Aq^j to all orders. Further numerical support for 
this interpretation of the quenched pseudoscalar double-hairpin comes from the shape of the 
double-hairpin propagator as a function of Euclidean time. As discussed below (see also Ref. 
||22||), this time dependence can be fit very well to the functional form (1 + m^t) exp{—mTjt) 
expected from a mass insertion vertex surrounded by two propagators of mass m^r. 

B. Methods 

The ability to study the double-hairpin diagrams relevant to the OZI rule has been greatly 
improved by two recent developments in lattice QCD methodology. The global source tech- 
nique (which we refer to as the "allsource" method) was introduced several years ago for the 
purpose of studying the rj' mass and the ?7a(1) anomaly ||2l|. In this method, the quark prop- 
agator is calculated from a sum of identical unit color-spin sources located at all space-time 
points on the lattice. If this allsource propagator is contracted over color indices at a given 
site, the result is a gauge invariant term corresponding to a closed quark loop originating 
from that site, plus a very large number of gauge-dependent open loops. The latter terms 
tend to cancel due to their random phases, allowing a determination of closed loop averages 
and loop-loop correlators (double-hairpins). The other recently developed technique which 
has greatly improved the accuracy of the results for both double-hairpin calculations and for 
other chiral studies with Wilson-Dirac fermions is the Modified Quenched Approximation 
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(MQA) p2[, which provides a practical resolution of the exceptional configuration problem 
that has long plagued such calculations. This method identifies the source of the exceptional 
configuration problem as the presence, in some gauge configurations, of exactly real eigen- 
modes which are displaced into the physical mass region by the artificial chiral symmetry 
breaking associated with the lattice Wilson-Dirac operator. By systematically identifying 
these real eigenmodes and calculating their contribution to the quark propagators, the cor- 
responding propagator poles can be extracted and moved to zero quark mass (where they 
belong). This MQA procedure has been applied to both the allsource propagators for double- 
hairpin calculations as well as to valence quark propagators. The resulting MQA-improved 
propagators have recently been used in an extensive study of quenched chiral logs and their 
relation to the rj mass and the anomaly [^. As a part of this study, the size and 

time-dependence of the pseudoscalar r( double-hairpin diagram was calculated, using the 
allsource method. Since the pole-shifting procedure has already been applied to the quark 
propagators, it requires very little additional effort to investigate the vector, axial-vector, 
and scalar double-hairpins which determine the spin-parity pattern of OZI mixing. The 
results we present here are from a set of 300 quenched gauge configurations on a 12^ x 24 
lattice at /? = 5.7. In the study of Refs. P^ , p5| , both naive Wilson and clover improved 
quark actions were studied. It was found that, at /3 = 5.7 with the Wilson action, sub- 
stantial lattice spacing effects suppressed the pseudoscalar double-hairpin, giving a smaller 
than expected value of A^^j =(0.27 GeV)^ for the double-hairpin contribution to the rl 
mass versus the value (0.49 GeV)^ extracted from weak-S'[/(3)-breaking mass formulas [jl4 



A much more satisfactory result is obtained from the clover improved quark action. With 
a clover coefficient = 1.57, the pseudoscalar double-hairpin gives Aq^j = (0.41 GeV)2. 
For the calculation of OZI-violating amplitudes, we will therefore use the clover improved 
quark action only; we also use the physical charmonium IS-IP splitting to set the scale 



1.18 GeV) for /3 = 5.7 when we quote lattice results in physical units p3 
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C. Results 



Using the method described in the previous Section, we have calculated the double- 
hairpin contribution to matrix elements of the form 

(g'(|/)rg'(y) q{x)rq{x)) (3) 

with Hermitian operators generated by the choices F* = ^75 (pseudoscalar), = 7^^, /i = 
1,2,3 (vector), F* = 7'^75, /i = 1,2,3 (axial vector) and F* = 1 (scalar) (the antisymmetric 
tensor a'^^ does not explore new states: it also has axial vector quantum numbers). As 
in standard hadron spectroscopy, we Fourier transform the space-time propagator over 3- 
dimensional time slices at zero 3-momentum and study its time-dependence. A particular 
advantage of the allsource method is that the Fourier transforms can be performed over both 
ends of the meson propagator, unlike the usual case of a fixed local source where only one 
end can be transformed. This provides an improvement in statistics which is quite important 
for the success of the method. For the scalar double-hairpin matrix element, the expectation 
value of a single scalar loop is nonzero, and so a constant proportional to (0|gg|0)^ must be 
subtracted from the above matrix element to get the true correlator. 



TABLE II. Quenched lattice OZI- violating amplitudes. 



nonet 


(MeV)2 


Ai'^j (MeV)* 


0-+ 


+(407 ± 11)2 


~ +290 


1— 


< (220)2 


< 30 


1++ 


< (380)2 


< 60 


0++ 


-(1350 ±90)2 


~ -520 



* For the conversion from Aq'^j extracted from the lattice via Eq. (^ to Aq'^j for comparison 
to the amplitudes quoted in Table | based on mass matrices, see Ref. | |14| |. 
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Even without any detailed analysis, the overall empirical OZI pattern of Table | is strik- 
ingly confirmed by the lattice results. This is easily seen from the size of the various double- 
hairpin correlators. In Figs. we have plotted the double-hairpin correlators for the 
pseudoscalar, vector, axial vector, and scalar sources. All plots have the same scale for 
comparison. The calculations have been done for 9 different choices of quark mass. The 
data shown in the figures are from one of the lightest quark masses, for which the pion 
mass is about 300 MeV {m.,^a = 0.266 ± 0.004). The results quoted in Table f| are chirally 
extrapolated to the physical pion mass. The errors in Figs. 6-9 and in Table f| are statis- 
tical only. By far the largest and longest-range correlator is the pseudoscalar correlator of 
Fig. ^. This is expected for two reasons: the anomaly introduces a large double-hairpin 
vertex responsible for the large r]' mass, and, as explained above, in the quenched approx- 
imation the external qq meson propagators on either side of the double-hairpin vertex are 
light Goldstone bosons. The results extracted from Fig. ^ have been reported in Ref . 



Compared to the very strong pseudoscalar double-hairpin, the vector and axial vector 
double-hairpins of Figs. |^ and ^ are dramatically suppressed, consistent with the empirical 
observations described in Section |I[ Since quenched lattice QCD gives reasonable values for 
the three-point functions associated with the meson virtual loop processes depicted in Fig. 
^, these results provide not only a first derivation of the OZI rule from QCD, but also a 
dramatic example of the evasion in QCD of the "second order paradox" described in Section 
I and a confirmation of the fact that in a complete calculation a conspiracy of the type 
described in Section |I|must occur. (Of course the results reported here include not only the 
meson loop contributions but also the other time orderings of the double-hairpin graphs of 
Fig. 0(b).) We in fact see no significant signals in the vector and axial vector channels and 
so report in Table y only one standard deviation upper bounds. 

As described in Section y and illustrated in Fig. ^ if the conspiratorial cancellation 
amongst meson loops is associated with ^Pq pair creation, one would expect Aq^^i be 
very large. Fig. |^ shows this behaviour: after taking into account the heavier mass of 
the scalar meson (about 1.3 in lattice units [2^), we find that the scalar OZI amplitude is 
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comparable in size to the pseudoscalar amplitude but of the opposite sign (see Table ||) . A 
full amplitude Aqzi in general has glueball, instantaneous, and loop contributions, and in a 
given amplitude, any or all of these components might be important. (Recall, for example, 
that while the loop contribution to Aq^j is believed to be small ||2^, the full Aq^j is large.) 
That the measured A'q^j is actually consistent in sign and magnitude with the hadronic 
loop contribution predicted by the quark model has interesting implications which we will 
discuss below. A large and negative Aq^^j has been previously reported in Ref. [p^] . 

To obtain the quantitative results for the OZI mixing amplitudes quoted in Table f|, 
we carried out an analysis similar to that used to obtain the rj' mass from the pseudoscalar 
double-hairpin p2| , |23[| . For that case, the time-dependence of the pseudoscalar double- 
hairpin correlator corresponding to Fig. ^(b) was found to be quite well described by a 
"double-pole" form consisting of a p^-independent double-hairpin insertion between a pair 
of meson propagators (see also Ref. [1^). In momentum space 

A.(p) = - fp^^A'ozj^^fp (4) 

pZ _|_ ^2 pZ _|_ 

where fp is the vacuum-to-one-particle matrix element 

fp = (0|g*7'gk(p)) (5) 
and Aqzi is the (mass)^ version of the Aqzi defined previously (called ml in Refs. 



2^ , |23[| ). This gives a time-dependent double-hairpin correlator at zero 3-momentum of the 



form 

A,(p = 0; t) = - ^^^1^(1 + m.t)e-™-* + (t ^ (Na - t)) (6) 
to be compared to the usual valence quark (e.g., isovector) correlator corresponding to Fig. 



0(a) 



A.(P) = fp-r^fp (7) 



which gives 
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A,(p = 0; t) = -iLe— + (t^ {Na - t)) 
2m^ 



(8) 



(The relative sign of Eqs. (^) and (|^) is tricky; with our convention a positive Aqzi makes 
a positive contribution to the (mass)^ of a state.) Since the values of fp and m^r can be 
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FIG. 6. The pseudoscalar double-hairpin correlator. Note the scale: this amplitude is negative. 
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FIG. 7. The vector double-hairpin correlator. 
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FIG. 8. The axial vector double-hairpin correlator. 



3 - 

c 
"Q. 




1 2 3 4 5 6 7 

time (iattice units) 



FIG. 9. The scalar double-hairpin correlator. 

separately determined from fitting Eq. to the valence quark correlator, the double- 
hairpin vertex insertion Aq/j can be determined by a one-parameter fit of (P) to the overall 
size of the double-hairpin correlator. A similar analysis of the scalar double-hairpin led to 
the result quoted in Table |I|, while for the other channels such analyses provided the quoted 
upper bounds for the very tiny mixing amplitudes in these channels. 
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IV. CONCLUSIONS 



The most straightforward conclusions of this work are that QCD can explain the OZI 
rule in channels where it is observed and that it predicts that Aq^^j is large and negative 
p6| . This supports the quark model's explanation of the dynamical suppression of the 
typical scale of hadron-loop-induced OZI violation below 1/Nc expectations, and in so doing 
provides further evidence for the standard ^Pq pair creation amplitude, since this is the 
critical feature which produces this result. 

While the precise consequences are unclear, the implications for phenomenology are 
serious. With Aq^^j large, the lightest scalar meson nonet (the IP states) will be close to 
the SU{3) limit. We may therefore expect an octet of scalar mesons in the 1400 MeV range 
with the other IP states while the nearly singlet scalar state will be substantially lower in 
mass. Thus the usual assumption of phenomenological analyses that this region will contain 
the unmixed nonet of IP states and the scalar glueball is incorrect. For example, this region 
might well contain the isoscalar state of the 2P nonet. In addition, since ^oz/ is comparable 
to the IP — 2P splitting, there is no reason to assume that either the IP or 2P singlet's 
properties can be related by nonet symmetry to those of its octet. The net effect is that 
the definitive extraction of the glueball state from the scalar meson spectrum may be quite 
difficult. 

Given the importance of this task, it is certainly worthwhile to study the scalar mesons 



more carefully in the light of this result . On the lattice it might be possible to obtain the 
matrix of OZI-violating amplitudes connecting the cj-like and 0-like IP and 2P states; in 
models the low-lying scalar meson spectrum can be studied including the effects of a strong 
annihilation channel. 

Perhaps most critical is to use quenched lattice calculations of the mixed propagators 
from quarkonia to glueballs [^] to help resolve the scalar meson OZI violation reported 
here into the contributions of qq'q'q intermediate states, purely gluonic intermediate states 
associated with "true double-hairpin" graphs, and instantaneous contributions. Ultimately, 
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quenched lattice calculations of three-point functions could directly check the predicted 
negative loop contributions to Aq^^j by measuring the vertex functions which are the "raw 
ingredients" of the quark model calculation. In particular, in other than the 0^+ channel, 
one should see the required magnitudes and opposite signs of the virtual P-wave decays to 
two i = mesons and the S'-wave decays to one i = 1 and one £ = meson required to build 
up the near cancellation that is at the heart of the quark model mechanism. In contrast, for 
0^^ mesons these channels should have the same sign. 

V. DISCUSSION 

The results described here clearly have serious implications for the spectroscopy of O^"*" 
states, and define the series of investigations described above required to clarify the physics 
behind A'q^j. Such investigations are not only important for their impact on phenomenology, 
however. They are also important because our results highlight other more fundamental 
questions raised long ago by Witten on the apparent conflict between the instanton 
solution of the t]' mass (i.e., [/^(l)) problem and the large Nc limit. The quark model 
mechanism for the loop contributions to A'q^j is based on large A^^^. While our discussion 
of the ^Pq model has focused on its prescription for the quantum numbers of the created 
qq pair, it is also an essential ingredient of the model that this pair creates {qq') + (q'q) 
and not (qq) + {q'q') mesons, i.e., that it respects the OZI rule at tree level. The physical 
picture behind this feature of the model is that pair creation (at order 1/Nc) occurs by the 
breaking of the color flux tube connecting q and q. More generally, as mentioned above, this 
limit provides the only known field-theoretic basis for the success of not only the valence 
quark model, but also of Regge phenomenology, the narrow resonance approximation, and 
many of the systematics of hadronic spectra and matrix elements p!3|j28| -p0| . In contrast, 
it is widely believed that the Ua{^) problem is solved through instanton contributions to 
the axial anomaly. However, as emphasized by Witten, instantons vanish like e~^'= and 
so do not appear in the large N^. expansion. "Insofar as [instantons play] a significant 
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role in the strong interactions, the large Nc expansion must be bad. It is necessary to 
choose between the two." Note that these arguments draw an important distinction 

between semiclassically calculated instanton effects, which vanish like e~^'=, and more general 
topological gauge fluctuations, which can contribute at order 1/Nc to m^/. The real issue 
is not whether there are large fluctuations of FF in the QCD vacuum, but whether these 
fluctuations arise as local semiclassical lumps with quantized winding numbers or simply as 
a result of the generically large gauge fluctuations of a confining vacuum. 



To place this conflict in context, recall Eqs. (1) and (2). From Section |T| it is apparent 
that the amplitude for any of Nj massless qq pairs to annihilate to any other pair is the 
same, i.e., that AT does indeed have the form of the Nf = 2 matrix shown in Eq. (2). 
As explained earlier, this is consistent with the 't Hooft instanton interaction since the 
"scattering" amplitude S in Eq. (1) contains a contribution —A from instantons. Thus to 
leading order in A the decomposition of Eqs. (1) and (2) is general and the analyses of OZI 

violation in Refs. including that in the pseudoscalar sector are valid. It follows 

that from a purely phenomenological perspective it is irrelevant whether or not there is an 
instanton contribution to hadronic physics: a phenomenology with Aq^j 7^ is "legal" in 
any case, since the anomaly allows a resolution of the Ua{1) problem with or without 



ri| , p!2[| instantons. What remains unclear is the physics behind the annihilation amplitudes. 
Since a lattice simulation sums over all paths, it contains the instantons as tunnelling events 
between classical vacua, but the Feynman diagrams of QCD, which represent the quantum 
corrections around these vacua, are incapable of representing instanton physics. Thus if 
instantons are important in QCD, Feynman diagrams would have to be supplemented by 



effective interactions (like the 't Hooft interaction). As noted by Witten [|r^, the foremost 



victim of the failure of Feynman diagrams implied if instantons are important would be 
the large A^^^ expansion, since it assumes that all-orders properties of the QCD Feynman 
diagrammatic expansion are properties of QCD. 

The observations reported in this one more item to a growing and 

closely linked set of issues where the physics of instantons and the physics of large Nc 

22 



confront each other. Assuming that confinement and the Nambu-Goldstone mechanism 
are properties of the all-orders Feynman diagrammatic expansion of QCD, the large Nc 
expansion provides a consistent framework embracing all strong interaction phenomena. 
Among these phenomena are the hadron spectrum for all flavors of hadrons (including the 
1/A^c-suppressed hadronic widths which seem to be critical to Aq^^j), the OZI rule (now 
including Aq^^j), and the qq condensate. As Witten argued long ago [ill], given the Ua_{1) 



anomaly and conflnement, the large Nc limit is also capable of explaining the t]' mass at 
order 1/Nc without instantons. 

While its limited range of applicability makes it somewhat less attractive for phenomenol- 
ogy (instantons offer a competing explanation only for the properties of the lightest SU{3) f 
hadrons), the instanton picture |^ has received strong support from recent lattice results 
]. Measurements of the topological charge |Tl| of "cooled" gauge conflgurations show 



that in such circumstances this charge is quantized and localized as expected for instantons. 
Moreover, the zero-modes of the Dirac operator associated with the solution of the Ua{1) 
problem and the near-zero-modes associated with the qq condensate are also localized and 
in "cooled" conflgurations can be associated with these same instantons. The lattice results 
on these and other hadronic properties are consistent with the instanton liquid model P2 



Since, as argued by Witten, conflnement can replace instantons as the source of the 
Ua{1) anomaly and since conflnement can also produce a space-time localization of the 
origin of the t]' mass and of the qq condensate, in our view the true origin of these effects 
remains unsettled. The results of this paper may help to resolve this situation since for Aq^^j 
the two competing pictures lead to mechanisms that are very distinct. Flux-tube-breaking 
pair creation, a prototypical large Nc phenomenon, led to the prediction that the hadron 
loop contribution to ^oz/ large and negative as found here. Moreover, as stated in the 
beginning of this paper, quark models, with their conflned constituent quarks, naturally 
generate a large positive Aq^j In this case the loop contribution should be typically 
small 1^ , and the large positive quark model amplitude is associated with an instantaneous 



interaction. Instantons, through the instantaneous 't Hooft interaction, would lead to a 
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superficially similar pattern of OZI violation: a large positive Aq^j and a large negative 
A^Qzi- However, the origins of the large negative A^q^i are very different in the two cases: 
the instanton A^^i is associated with an instantaneous contribution while the quantitative 
similarity between the quark model prediction and our measured ^qzt suggests that this 
amplitude is associated instead with the meson loop contributions. 

Our result thus favors the large Nc and not the instanton interpretation of the solution 
to the 77' mass problem. Nevertheless, while suggestive, the quark model prediction is not of 
sufficient quantitative accuracy for this conclusion to be reliable. Fortunately, with recent 
advances in lattice methods and in computing power, we believe that the results we have 
described here can not only be improved but also understood more deeply. In particular, 
through the program we described of decomposing the OZI-violating amplitudes into their 
component parts, it should be possible to define the mechanism driving Aq^^j. We also 
believe it will be particularly fruitful to define and test confinement-based interpretations 
of the lattice results on such quantities as the topological susceptibility, the localization of 
zero modes, the correlation function of the topological charge operator, and the space-time 
association of the qq condensate with the topological charge. Through such studies, the 
confiict between large Nc and instanton physics can at last be resolved. 



24 



ACKNOWLEDGEMENTS 

We are grateful to Stephen Sharpe and Thomas Schaefer and to Chris Michael for alert- 
ing us to a serious sign error in the first version of this paper and to important references 
which had escaped our attention. This work was supported by DOE contract DE-AC05- 
84ER40150 under which the Southeastern Universities Research Association (SURA) oper- 
ates the Thomas Jefferson National Accelerator Facility. The work of H.B. Thacker was 
supported in part by the Department of Energy under grant DE-FG02-97ER41027. 

REFERENCES 



[1] G. Zweig, CERN Report No. 8419 TH 412, 1964 (unpublished); reprinted in Developments in 
the Quark Theory of Hadrons, edited by D. B. Lichtenberg and S. P. Rosen (Hadronic Press, 
Massachusetts, 1980). The history of the discovery of the quark model (or "aces") as seen 
by Zweig is related in "Baryon 1980", Proceedings of the IVth International Conference on 
Baryon Resonances, ed. N. Isgur (Toronto, 1980), p. 439. 

[2] S. Okubo, Phys. Lett. 5, 165 (1963); Phys. Rev. D 16, 2336 (1977); J. lizuka, K. Okada, and 
O. Shito, Prog. Th. Phys. 35, 1061 (1966); J. lizuka, Prog. Th. Phys. Suppl. 37-38, 21 (1966). 

[3] A. De Rujula, H. Georgi, and S. L. Glashow, Phys. Rev. D12, 147 (1975). 

[4] N. Isgur, Phys. Rev. D12, 3770 (1975); D13, 122 (1976). 

[5] H. Fritzsch and P. Minkowski, Nuovo Cimento 30A, 393 (1975). 

[6] There is an interesting subtlety associated with this decomposition. Fig. ^(b) includes time- 
orderings in which the pair creation occurs before the annihilation (see Fig. |^) corresponding 
to meson loop processes. However, the complete set of meson loop processes produce not only 
Fig. ^(b) but also Z-graphs of Fig. ^(a) and graphs with closed qq loops inserted into Fig. §(a). 
Thus a consistent treatment actually requires that S be not just the simple valence graph 



25 



shown, but also include graphs with internal qq loops Q. 

[7] M. Gell-Mann, R.J. Oakes, and B. Renner, Phys. Rev. 175, 2195 (1968); S. Weinberg, Phys. 
Rev. D12, 3583 (1975). 

[8] Y. Nambu, Phys. Rev. Lett. 4, 380 (1960); J. Goldstone, Nuovo Cim. 19, 154 (1961). 

[9] A.M. Polyakov, Phys. Lett. 59B, 82 (1975); Nucl. Phys. B121, 429 (1977); A.A. Belavin, A.M. 
Polyakov, A. Schwartz, and Y. Tyupkin, Phys. Lett. 59B, 85 (1975); C. Callan, R. Dashen, 
and D.J. Gross, Phys. Lett. 63B, 334 (1976); R. Jackiw and C. Rebbi, Phys. Rev. Lett. 37, 
172 (1976). 

[10] G. 't Hooft, Phys. Rev. Lett. 37, 8 (1976); Phys. Rev. D14, 3432 (1976). 
[11] E. Witten, Nucl. Phys. B149, 285 (1979); 256, 269 (1979). 
[12] G. Veneziano, Nucl. Phys. B159, 213 (1979). 

[13] G. 't Hooft, Nucl. Phys. B72, 461 (1974); E. Witten, Nucl. Phys. B160, 57 (1979). 

[14] The extraction of the OZI amplitudes from the data can be complex. If Aqzi is small, uj — (p- 
like mixing is small as is mixing to excited nonets and glueballs, so Aqzi may readily be 
extracted from the relation mi=o = mi=i + 2Aozi- Note that since this extraction is in the 
SU (2) f sector, issues of 5'[/(3)-breaking do not arise. Also note that the relation between 
Aqzi and Aqzi defined by the analogous (mass)^ formula rnj^Q = rnj^^ + 2Aozi is simply 
Aqzi = Aozi/2m[=i at the light quark mass scale. When Aqzi is large, mixing and SU{3) 
breaking can become important. In such circumstances, no simple comparison to the data can 
be made without further assumptions. Since excited nonets will normally be 0.5 to 1.0 GeV 
away, when Aqzi is as large as it is in the pseudoscalars and scalars, the validity of ignoring 
internonet mixing is certainly questionable. If one nevertheless makes this assumption, then in 
the S'C/(3) limit one would have mo = 'mi=i+3Aozi or alternatively ttt-q = rnj^^+SAozi where 
niQ is the mass of the SU (3) singlet meson, leading to the relation Aqzi = Aozi/['mo + ^i=i]- 
For weak SU{3) breaking these formulas simply become mg = 2mj^i — m'j^Q + SAqzi and 

26 



Q = 2m^_ 1 — ni'f^Q + 3 Aqzi where now mg and rn'j^^ are the masses of the mainly singlet 
and mainly octet mesons (e.g., the rj' and the rj, respectively). It is this latter formula that is 
commonly used in defining A.q^j in the Witten-Veneziano formula [^.12] and which leads to 
SAl^zi =(0.85 GeVf. These broken SU{3) relations lead to Aqzi 



A 



ozi- 



The inaccuracy of these formulas from both internonet mixing and higher order corrections 
in SU (3)-breaking are not obviously small. Thus, given the dynamical assumptions required 
to relate a large value of Aqzi to observed masses, the empirical value of Aq^j (or Aq^j) is 
quite uncertain (as is A^^j and A^^^) and therefore only semiquantitative statements about 
these amplitudes can be made at this time. To convert the lattice propagator data for Aqzi 
to Aqzi in Table || for comparison to the Aqzi quoted in Table |, we note that in the single 
flavor case considered, if we once again make the assumption of small mixing with other states, 
the mass and (mass)^ formulas are m = mj=i + Aqzi and = rnj^i + Aqzi- (Note that 
we have identified the mass without the annihilation amplitude as being the 1 = 1 mass that 



would be found for two or more flavors.) It follows that Aqzi = y"^?=i + -^oz/ — "i/=i. 
For small Aqzi this relation gives Aqzi — Aqzi /'^^i=i and (as in our extraction of these 
quantities from the data) in this case our assumption of small mixing is justifled. For large Aqzi 
such an assumption is probably not good; however, the lattice relation can be systematically 
improved by measuring additional elements of the mass and (mass)^ matrices should such an 
improvement be justifled by an improvement in the relevant experimental data. 

[15] H.J. Lipkin, Nucl. Phys. B291, 720 (1987); Phys. Lett. B179, 278 (1986); Nucl. Phys. B244, 
147 (1984); Phys. Lett. B124, 509 (1983). 

[16] There is a loophole in this argument. It is possible that instead there could be a conspiracy 
between the hadronic loop processes and other sources of OZI-violation. 

[17] P. Geiger and N. Isgur, Phys. Rev. D 44, 799 (1991); Phys. Rev. Lett. 67, 1066 (1991); Phys. 
Rev. D 47, 5050 (1993); P. Geiger, ibid. 49, 6003 (1993). 

[18] L. Micu, Nucl. Phys. BIO, 521 (1969); A. Le Yaouanc, L. Oliver, O. Pene, and J.-C. Raynal, 



27 



Phys. Rev. D 8, 2233 (1973); Phys. Lett. B71 , 397 (1977); ibid. B72 , 57 (1977); W. Roberts 
and B. Silvestre-Brac, Few Body Syst. 11, 171 (1992); P. Geiger and E.S. Swanson, Phys. Rev. 
D 50, 6855 (1994); R. Kokoski and N. Isgur, Phys. Rev. D 35, 907 (1987); Fl. Stancu and 
P. Stassart, Phys. Rev. D 38, 233 (1988); 39, 343 (1989); 41, 916 (1990); 42, 1521 (1990); 
S. Capstick and W. Roberts, Phys. Rev. D 47, 1994 (1993); Phys. Rev. D 49, 4570 (1994); 
E.S. Ackleh, T. Barnes, and E.S. Swanson, Phys. Rev. D 54, 6811 (1996); and for a recent 
and detailed description of the implications of the ^Pq model for mesons, see T. Barnes, F.E. 
Close, P.R. Page, and E.S. Swanson, Phys. Rev. D 55, 4157 (1997). 

[19] C. Schmid, D.M. Webber, and C. Sorensen, Nucl. Phys. Bill, 317 (1976); E.L. Berger and 
C. Sorensen, Phys. Lett. 62B, 303 (1976). 

[20] Unfortunately, Ref. |17| did not check the loop contribution to ^oz/ since the Ua{^) anomaly 
guaranteed that OZI violation in this sector would be large so that there was no second order 
paradox to evade. However, given the identified mechanism of the suppression of the loop 
diagrams, there is every reason to suppose that the loop contribution to Aq^j is as small as 
in any other nonet with J^'-' ^ 0^^. 

[21] Y. Kuramashi, M. Fukugita, H. Mino, M. Okawa, and A. Ukawa, Phys. Rev. Lett. 72, 3448 
(1994). 

[22] W. Bardeen, A. Duncan, E. Eichten, and H. Thacker, Nucl. Phys. B (Proc. Suppl.) 73, 243 

(1999) . 

[23] W. Bardeen, A. Duncan, E. Eichten, and H. Thacker, Nucl. Phys B (Proc. Suppl.) 83-84, 215 

(2000) . 

[24] A.X. El-Khadra, G.M. Hockney, A.S. Kronfeld, and P.B. Mackenzie, Phys. Rev. Lett. 69, 729 
(1992). 

[25] We have determined the p, oq, and a\ masses from chirally extrapolated fits to our data 
on the valence (i.e., isovector) propagators associated with Fig. |2|(a), which give ~ 0.8 



28 



GeV, rriai ~ 1.4 GeV, and mag — 1.6 GeV. The latter two masses seem to be too high. It 
is quite possible that they are especially sensitive to quenching since both the ai and the uq 
have (or are predicted to have) very strong S-wave couplings to nearby thresholds. (See, for 
example, R. Kokoski and N. Isgur, Phys. Rev. D 35, 907 (1987); N. Isgur, Phys. Rev. D 57, 
4041 (1998).) We have found no other recent lattice measurements of these masses close to 
the chiral limit (see, however, Ref. [^6|), or of other spectroscopic properties of the P-wave 
mesons, and consequently hope to report on their spectroscopy in a forthcoming publication. 
In any event, at the level of accuracy relevant to the results reported here, these 10% effects 
are not significant. 

[26] C. Michael, M. S. Foster and C. McNeile (UKQCD Collaboration), Nucl. Phys. B (Proc. 
Suppl.) 83-84, 185 (2000). 

[27] W. Lee and D. Weingarten, Nucl. Phys. B (Proc. Suppl.) 63, 194 (1998); Nucl. Phys. B (Proc. 
Suppl.) 73, 249 (1999); Phys. Rev. D 61:014015 (2000). See also Ref. |2|] on ^Pp-glueball 
mixing. 

[28] The Dual Topological Expansion is very closely related to the large Nc expansion. See G.F. 
Chew and C. Rosenzweig, Nucl. Phys. B104, 290 (1976); G. Veneziano, Proc. 9*^ Ecole d'Ete 
de Physique des Particules, Gif-sur-Yvette, 1977 (1978), vol. 2, p. 23. 

[29] R. Dashen and A. Manohar, Phys. Lett. B315, 425 (1993); B315, 438 (1993); R. Dashen, E. 
Jenkins and A. Manohar, Phys. Rev. D49, 4713 (1994); D51, 2489 (1995). 

[30] Elizabeth Jenkins and Richard F. Lebed, Phys. Rev. D 52, 282 (1995). 

[31] E. Witten, Nucl. Phys. B149, 285 (1979) on page 286. 

[32] For a review, see T. Schafer and E. Shuryak, Rev. Mod. Phys. 70, 323 (1998). 

[33] For a review, see J. Negele, Nucl. Phys. B (Proc. Suppl.) 73, 92 (1999). 

[34] John F. Donoghue and Harald Gomm, Phys. Rev. D 28, 2800 (1983). 



29 



